INVESTIGATIONS OF THE MIYATA SYNTHESIS 
TECHNIQUE 



by 



Richard Eric Carlson 



LIBRARY 

FAVAL POSTGRADUATE SCHOOL 
MOHTEREY, CALIF. 93940 



United States 
Naval Postgraduate School 




THESIS 



INVESTIGATIONS OF THE MIYATA SYNTHESIS TECHNIQUE 



by 



Richard Eric Carlson 



June 1970 



Jkit, document hcu> been approved fioA. pubtic Ae- 
Ze.cn , e and bate.’, Z£t> di^VvibuutLon Z& unZUjruXe.d . 




Investigations of the Miyata Synthesis Technique 



by 



Richard Er ic^Car Ison 
Lieutenant, United States Navy 
BSEE, Purdue University, 1966 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1970 



- 




h 



1 













CWGEADUATE SCHOOTJ 
Y < CALIF. 93940 



ABSTRACT 

These investigations generalize Miyata's synthesis of passive 
driving-point impedance functions by developing a step-by-step com- 
putational technique for augmenting a general driving-point impedance 
to insure that the real part of the augmented impedance is positive 
term by term. This goal has been accomplished and programmed under 
the provision that the real part of the original impedance is minimum 
reactive, i.e., has no zeros on the jo)-axis. The latter requirement 
is not restrictive since zeros on the jcu-axis can be removed a priori 
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I. INTRODUCTION 



7 8 

The Miyata ’ method of transformerless synthesis provides a method 
for separating a given impedance into additive components which indi- 
vidually satisfy certain special conditions that permit their realization 
in simple terms. This method of separation is based upon the even part 
of the given impedance and can be stated briefly as follows: Given a 

minimum reactive driving point impedance Z(s), find R(co) the real part 
of Z(s). Then split R(co) into a sum of even functions R^(cn), (p=0 , 1 , . . . ,n) 
each satisfying the requirements of ReZ^(jco). Find Z^(s) corresponding to 
each Rp(o)). The realization of Z(s) reduces to the realization of each of 
the Zp( s ) an d the series connection of the resulting one ports completes 
the synthesis. 

An inherent difficulty in the above technique lies in the require- 
ment that the Rp(to) t> e positive for all cn. To overcome this difficulty 
Miyata suggested the use of surplus factors. An example will illustrate 
this method. 

Let Z(s) be given as 



Z(s) 



2 

6s +9s+9 
2 

s +3s+4 



( 1 . 1 ) 



from which we obtain 



R(a>) 



6od^ - 6a)2-f 3 6 

B(co 2 ) 



= R 2 (oo) + + R Q (to) 



( 1 . 2 ) 



This so called M 
R^(cjd) is negative, 
augmented impedance , 



single-n split" method cannot be applied because 
Therefore use a surplus factor (s+1) to obtain an 



Z a (s). 



Z (s) 
a x ' 



s+ 1 
s+1 



Z (s ) 



6s^+15s 2 +18s+9 

s^+4 s 2 +7 s+4 



(1.3) 
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for which the numerator of ReZ (ja>) contains only positive terms. 

a 



R(cd) 



6cu^+30(D^+36 

B' (cjd 2 ) 



( 1 . 4 ) 



Obviously the use of surplus factors complicates matters since 

higher-order factors are often required. 

1 8 

At this point most authors 5 considering the technique, including 
Miyata, generally abandon the use of surplus factors without discussing 
any general method for their selection and proceed to describe various 
groupings of the individual positive and negative R^(cjd) which remain 
positive for all co. These grouping methods generally require factori- 
zation of the numerator with its attendant difficulties and no general 
method is available in the literature nor is it clear that such groupings 
are always possible. 

The purpose of this thesis is to investigate in detail the use of 
surplus factors to develop augmented impedances that are positive term 
by term. 

Chapter two is a general consideration of the Miyata method which 
describes the reasoning behind the technique and the mechanics of 
developing the electric circuit and component values from the original 
driving-point impedance when all terms in the numerator of ReZ(jco) are 
initially positive. The significant results of this section include: 

The circuit realization is transformerless, an important advantage 
over the Bott-Duffin technique. 

The complete circuit development consists of a simple unbalanced 
ladder network with series L's and shunt C's terminated in a resistance. 

By term grouping without factorization, the number of circuit elements 
required is less than for the modified Bott-Duffin for all driving-point 
(dp) impedance functions of second order or greater. The saving in 



e* 
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number of elements becomes quite dramatic for higher-order polynomials 
so that even though the order of the original dp impedance must often 
be increased by the use of surplus factors, large savings are available. 

In general, chapter two is included to provide sufficient incentive 
to undertake the task of developing a general method of selecting 
surplus factors of minimum order. 

In chapter three, starting with a dp impedance that is the ratio of 

two general polynomials of the same order in s and a general augmenting 

polynomial (surplus factor), the real part of Z(s) is developed. Since 

2 n 

the resultant coefficients of od must be non-negative, these coefficients 

are the required constraints that must be satisfied to insure that ReZ ( jo)) 

is negative term by term. These coefficient constraints constitute a 

set of non-linear inequalities in the variables 6^ (i=0 , 1 , . . . ,n) which 

are the coefficients of the augmenting polynomial (i.e., ^s+ . . .+b^s n ) . 

Also in this section, the non-linear inequalities developed are linearized 

9 2 

using a set of congruent transformations and solved by the simplex method. 

The difficulties involved in the inverse transformation are circumvented by 

considering the form of the solved constraints which are the coefficients 

of the real part of the augmented dp impedance. The process is completed 

by demonstrating how a knowledge of the form of the solved constraints 

allows one to find the numerator of the real part by simple long division. 

Determination of the denominator is straightforward and the problem is 

solved. 

A digital computer program for part of the routine described above is 
developed in chapter four, and a sample problem worked in section five. 

In summary, a method has been developed which will invariable find 
an optimal augmenting polynomial such that the resultant augmented dp 



9 



impedance has a real part which is positive term by term. The method 
is optimal in the sense that as many coefficients as possible of the 
real part are zero, thereby reducing the required number of circuit 
elements to a minimum. 

The value of this method lies in the fact that it will invariably 
produce a positive real part for the augmented impedance in a step-by- 
step fashion. In fact, the entire process could be programmed such that 
for any input minimum-reactive dp impedance the output would be the 
required circuit with accompanying element values. 

The single major disadvantage is that the augmented dp impedance may 
have a large number of terms which would consequently require a large 
number of circuit elements to realize (the circumstances that determine 
this possibility are discussed in chapter two). 
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II. GENERAL CONSIDERATION OF THE MIYATA SYNTHESIS TECHNIQUE 



The Miyata method for realization of a passive, minimum-reactive, 

1,7 



driving-point impedance proceeds as follows: 
no jcn-axis zeros or poles written 



For an impedance having 



Z (s) = 



+ n^ 
m 2 + n 2 



( 2 . 1 ) 



where m denotes an even polynomial in s and n denotes an odd polynomial 
in s . 

Its real part for s=jo) is 

" n 1 ^A 

ReZ(jo)) = 2 2 (2.2) 

or 

a -la 2 . 2n . .2 a 2n 

A +A,a) +. . ,+A cd A A cd A cd 

ReZ(jco) = 5 oT = J~ + o - + + y~ (2.3) 

B +B,cjd +. . .+B a) B(cd ) B(cjd ) B(cjd ) 

o i n 

Each of the terms on the right-hand side, when the coefficients 
^o , ^l , **’ , ^n are P os: *- t: ^ ve 9 ma y be shown to generate an impedance that 
can be realized without mutual inductive coupling. The series con- 
nections of such separate realizations yields the desired synthesis of 
Z (s) . 

To see why this is so, consider the component impedances of the real 
part. The first term has n double-order zeros at s= infinity; the second 
has one double-order zero at s=0 and n-1 zeros at s= infinity; the third 
term has two double-order zeros at s=0 and n-2 zeros at s= infinity; 
finally, the last term has all its n double-order zeros at s=0. There- 
fore, all of these terms have their double-order zeros located exclusively 
at s=0 and s= infinity. 
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When the real part of an impedance is zero at some joo-axis point, 
it need not necessarily follow that the entire impedance is zero. 

However, the jco-axis points under consideration are restricted to s=0 
and s= infinity and, since a minimum-reactive function must be purely 
real at these points, it follows that when the real part is zero the 
whole impedance must necessarily be zero. 

Next we use the Darlington^ theory to show that these impedances are 
completely developable in simple terms. The Darlington theory shows that 
in a cascade development of an impedance each component lossless network 
produces one set of zeros out of the total complement of zeros which the 
even part of that impedance contains. Therefore, the network generated 
by the first cycle provides one set of zeros and leaves a remainder whose 
even part contains the remaining even-part zeros. The next network in 
the chain provides another set and leaves a second remainder with the 
leftover even-part zeros - and so forth. 

Applying this reasoning to the real part of the first term of the 
expansion of the real past of Z(s), we can predict that the corres- 
ponding impedances will have a zero at s= infinity and its reciprocal 
a pole at this point. Removal of this pole as a shunt capacitor is one 
Darlington cycle. The admittance remainder after this pole at s= infinity 
is removed has an even part that still has n-1 double-order zeros at 
s= infinity since so far only one has been used up. This admittance 
remainder must have a zero at s= infinity and its reciprocal must have 
a pole there. Removal of this pole as a series inductor constitutes the 
second cycle. The remainder now has an even part with n-2 double-order 
zeros at s= infinity and is zero at s= infinity. Its reciprocal must 
have a pole at this point and its removal as a shunt capacitance completes 
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the third cycle. This process continues yielding series inductive and 
shunt capacitive branches of total number n, the number of double-order 
zeros possessed by the original real part. When the last of these zeros 
has been accounted for, the remainder is reduced to a constant because 
its real part has no more zeros. The complete development consists of a 
simple unbalanced ladder network with series inductances and shunt cap- 
acitances terminated in a resistance. Subsequent terms are similarly 
deve loped . 

Using this single-n split technique, the total number of reactive 
2 

elements is n , where n equals the highest degree of either the numerator 
or denominator polynomial in s. Comparison with the number of elements 
required for a Bott-Duffin synthesis is displayed in Table I along with 
the reductions due to the 1/2-n split technique which is subsequently 
discussed. The geometric character of the Bott-Duffin method compared 
with the arithmetic character of the Miyata method is strongly apparent 
in this table. Even though special devices are required to salvage the 
Miyata method when all the terms of the real part numerator are not 
initially positive, there appears to be considerable advantage to 
exploring the Miyata procedure where more elaborate impedance functions 
are involved. Before one gets too enthusiastic about the Miyata method, 
it would be well to consider the requirement that all the coefficients 
in the polynomial formed by the numerator of the real part of the 
driving-point impedance must be positive. 

A sufficient condition that the coefficients be positive can be 
obtained by recalling that a Hurwitz polynomial has positive coefficients. 
Thus if 



A(<x> 2 ) = A + ... + A <x> 2n 
0 n 



( 2 . 4 ) 
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has no zeros in the r ight-ha If oo plane, then A , . . . ,A are positive. 

U n 

Or if 



A(-s ) = - n 1 n 2 



(2.5) 



has no zeros in the left-half s plane, then the A^,...,A^ are positive. 

2 

Since the left-half s plane maps into a region of the s plane con- 
taining all points not more than plus or minus 45° from the joo-axis, we 
can say the coefficients in eq. 2.4 are positive so long as there are no 
even-part zeros of Z(s) closer than 45° to the joo-axis. 

Although the coefficients in 2.4 may still be positive when some 
even-part zeros lie closer to the joo-axis than this sufficiency con- 
dition permits, the chances of this desired condition being fulfilled 
when even-part zeros are relatively close to the joo-axis become remote 
and vanish completely when a single pair of even-part zeros lie upon the 
joo-axis . 

In a manner of speaking, we can say that a polynomial that has some 
zeros above the 45° line in the s plane may still yield all positive 
coefficients if there are a sufficient number of zeros below the 45° 
lines to produce a compensating effect. Therefore, in a situation in 
which the even part of Z(s) has some negative coefficients, we can augment 
the real part with appropriately chosen factors to bring about the desired 
compensation. 

The method of compensation fails if the even part of Z(s) has any 
joo-axis zeros, since the number of necessary compensating factors tends 
toward infinity as one or more pairs of even-part zeros approach the 
joo-axis. With such a compensating scheme the number of terms in the real 
part and hence the number of necessary elements increases rapidly. 
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One method to reduce the number of elements is considered now, the 

1 8 

1/2-n split. 5 Miyata points out that the detailed procedure for 
carrying out the basic method of decomposition may be varied by separating 
the real-part numerator polynomial into groups of terms rather than 
separate terms. 

Consider 



/a , A 2 IA 2 (n-k) v / 2 (n-k+ 1 ) 2 n. 

(A +A_o) + . . ,+A . a) ) (A a) + . . .+A a) ) 

.. .w. v o 1 n-k , n-k +1 n N 

ReZ(ja)) = 5 + 5 (2.6) 



B(c» 2 ) 



B(o> 2 ) 

where k=n /2 (n even) k=(n+l )/2 (n odd). 

S. N. Hunt , who refers to this decomposition as the 1/2-n split, finds 
that it yields a smaller total number of elements than any other division 
of the numerator polynomial into parts. It requires fewer elements than 
the Bott-Duffin method for any degree n equal to or greater than 2. The 
impedance generated by the first term is developable into a simple unbalance 
ladder with series L T s and shunt C f s (total number k) terminated in a 
remainder whose even-part zeros are the same finite non-zero ones contained 
in the even part of this impedance. The second term is developable into an 
unbalanced ladder with series C f s and shunt L f s (total number n-k+1) and a 
terminal impedance whose even-part zeros are those that are contained in 
the even part of the second impedance. To the terminal impedances which 
appear as remainder functions we can apply a similar treatment and continue 
until the ultimate remainder functions are resistances. 

The process of augmentation by surplus factors is also necessary in 

order that Miyata f s elegant method of determining an impedance from its 

7 8 

real part may be used. The method proceeds as follows: * 

m + n 

Z(s) = ~ t " 
m 2 + n 2 



(2.7) 
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and 



EvZ(s) = — 2 



m i m 2 ' n l n 2 



m n 



N M 

D(s) 



( 2 . 8 ) 



“2 2 

Next consider an auxiliary function 



Z'(s) = 



NCm^ + n 1 ') 
m 2 + n 2 



(2.9) 



Its even part is 



EvZ * (s ) = 



N (m ^ f “ n^'n^) 



( 2 . 10 



m n 



2 2 

Set the two even parts to be equal; 

EvZ 1 (s) = EvZ (s) 
which can be done provided 

m, f rru - n *n^ = 1 . 



( 2 . 11 ) 



( 2 . 12 ) 



1 2 1 2 

where the degree of m^+n^ 1 is that of m^+n^. Polynomials and 

n^ 1 can always be found to satisfy this requirements. Z'(s) however will 
have a numerator of degree two or more higher than its denominator and is 
therefore not positive real (p.r.). To find the p.r. function Z(s), carry 
out a long division on Z'(s) until the degree of the numerator does not 
exceed the degree of the denominator. That is. 



N(m_ 1 4- n ’ ) m + n 

Z' (s) 1 — - q(s) + 

rru + n^ rru + n rt 



(2.13) 



where q(s) is the quotient polynomial and m^+n^ is the remainder 



Let Z (s) be 
r 



m + n 

Z ( S ) = r r 
r m_ + n» 



(2.14) 



so that 



Z(s) = Z r (s) + q(s) 



(2.15) 
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Recalling that EvZ f (s) = EvZ(s), 

EvZ(s) = EvZ r (s) 4 Evq(s). (2.16) 

If the polynomial q(s) contains an even power of s then the right 
hand side of eq. 2.16 tends to infinity as s=jo) tends to infinity. On 
the other hand, EvZ(s) g _j^ = R(cq) is bounded at infinity. Therefore q(s) 
must be an odd polynomial and as a result 

EvZ^(s) = EvZ(s). (2.17) 

Since the numerator of Z^(s) does not exceed the denominator Z^(s) 
has no pole at infinity or at any other finite frequency on the jo)-axis 
since its denominator is that of Z(s). 

Briefly Z^(s) has no poles in the right-half plane or on the jco-axis, 
and its even part for s=jo) is always non-negative. Therefore Z^(s) is 
positive real. To illustrate: 



Given ReZ^Qo)) = 



1 4 o> 



(2.17a) 



„ , v _ 2/3 s + 4/3s + 1 

Z o (s) “ 2 

s 4 2s 4 2s 4 1 



(2.17b) 



from any method whatsoever. Suppose we want an impedance with the 

2 6 2 

real part ReZ^(jo)) = cd / lH-co which is 2.17a multiplied by cd and implies 

2 

m=-s . 

2 

Then F(s) equals eq. 2.17b multiplied by -s and the long division 
takes the form 



-2/3s 



3 2 

s + 2s 4* s 4 1 



-2/3s£ - 4/3s^ - s 2 
-2/3s - 4/3s - 4/3s 



l/3s 4 2/3s 



and 



„ , X _ l/3s 2 4 2/3s 

z i (s) " T~7 2 ~~ 



s + 2s + 2s + 1 
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2 n 6 

Using this procedure any Z^Cs) = cu /(I + o> ) n=l,2,... can be 

found from Z q (s). If we want an impedance with the real part 

2 4 6 

ReZ(jco) = A q + AjCO + A^cd + A^cn 
Z(s) may be generated as 



A q Z o (s) + A 1 Z 1 (s) + A 2 Z 2 (s) + A 3 Z 3 (s) 



where continuing as above 



„ , x 2/3s^ + 1/3 s 

Z 2 (s) “ 1 2 

s + 2 s + 2 s +1 



, t s _ s 3 + 4/3s 2 + 2/3s 
"o (s; o 

s + 2 s + 2 s + 1 



s f° un< ^ by multiplying Z^(s) by -s and Z^(s) found by 

2 

multiplying ^ 2 ( 3 ) by -s and performing the indicated long division. 
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Ill . GENERAL STEP-BY-STEP PROCEDURE FOR THE SELECTION OF SURPLUS FACTORS 



The topic of this chapter is the step-by-step procedure for develop- 
ing the real part of an augmented impedance that is positive term by term 

from a general n^-order, passive, minimum-reactive, driving-point 

t h 

impedance augmented by a general m -order augmenting polynomial. 



A. CONSTRAINT DEVELOPMENT 

This section is devoted to the development of the required constraints 
to insure a positive real part for the dp impedance and begins with a 
general impedance function and augmenting polynomial as follows: 



Consider first the general driving point impedance Z(s) 



N 

E 



Z (s) 



= q=0 a q S = N(s ) 
- N q - D(s) 

2 p q 

q=0 



( 3 . 1 ) 



The general augmenting polynomial P(s) 



PCs) = 2 6 s 

r=0 



( 3 . 2 ) 



The product of Z(s) and P(s)/P(s) is the augmented driving point 



impedance Z (s) ; 

a 



N T 

£ (a s q £ 6 s r ) 

z (s) - 3-° r :° - 

a v / N T 

2 (3_ s<1 2 6 s r ) 

q=0 q r=0 



Z a< S > = N 



N 

T 

2 a 6s 
q=0 qr 
r=0 



q+r 



m ^ 1 + n ^ 1 
m 2 ’ + n 2 7 



E 3 a 6 r sq+r 

q=0 q r 
r=0 



( 3 . 3 ) 



( 3 . 4 ) 
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The even part of the augmented driving point impedance is 



i „ t _ i _ i 

^ 1 ^ I X7l/\ 

„ _ , N 12 12 N'(s) 

EVZa<s) , 



(3.5) 



where 



(W-T)/2 

T 

m ' = E a 0 6 s 2u 
1 _ n 2u-r r 

u=0 

r=0 



(N4-T- 1)/2 
T 

n. ' = E a 0 , , 6 s 

1 2u+l-r r 

u=0 

r=0 



2u+l 



(NfT)/2 
T 

m ' = E 3 0 6 s 

2 . 2u-r r 

u=0 

r=0 



(W-T-l)/2 

T 

n ' = E 3 0 , t 6 s 
2 n K 2u+l-r r 

u=0 

r=0 



2u 



(3.6) 

2u+l 



and noninteger values of (N4-T)/2 and (N+T-l)/2 are rounded to the 
next lowest integer. 



Expanding the even part of Z (s) 

a 



(M-T)/2 

T 



2u 



(N4-T)/2 

T 



EvZ (s) = Za 0 6s Z 3 0 6s 

a . 2u-r r _2u-r r 

u=0 u=0 



2u 



(N+T-l)/2 

T 



2u+l 



(NT- l)/2 
T 



r=0 



r=0 



Z a 0 , , 6 s Z 3o j-i 6 s 

_ 2u+l-r r . K 2u+l-r r 

u=0 u=0 



2u+l 



r=0 



r=0 



D'(s) 



(3.7) 



Combining products , 



(NKD/2 (NhT- 1)/2 (3.8) 

EvZ (s) = Z s 2(u+v) a ? 3 ? 6 6 - Z s 2 ( u+v+1 ) a g 6 6 

3 u ,v=0 2v ' r 2u_q r q u ,v=0 2vfl-r p 2u+l-q r q 



q,r=0 



q> r= o 



D'(s) 



Combining the difference, 



(N4-T)/2 

Evz a ( s ) = E s 2( ^\s <c^ B 2 -s 2 ^ B 2 ) 



(3.9) 
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Since ReZ (jo)) = EvZ (s) and the Miyata method requires that 

a as ~ * jco 

the real part be positive term by term, constraints must be developed 
to insure this requirement. First, consider the denominator D*(s) above 

D ' (s) = m 2 ' 2 - n 2 ' 2 = D(ja))D(-ja>) = [ D ( joo) | 2 (3.10) 



which is never negative and is therefore of no concern. 

Thus the constraint is reduced to 

N* (a) 2 ) = m. f m ' - n- *n ' . =0 for 0=^jd=<» 

12 12 s=jo) 



(3.11) 



The coefficient of each term in eq. 3.9 must be non-negative which 

requires that 

(N+T-l)/2 

(NfT)/2 

T 

S s 4mf2 6 6 (a, 0 3, - a. .. 9 3 9 , . ) = 0 (3.12) 

_ r q 4nd-2-2u-r 2u-q 4m+l-2u-r 2u+l-q 

m=0 

u=0 

r,q=0 



and 



(N+T-2)/2 

(NfT)/2 

T 

„ 4mH-4 , . , * ^ „ 

Z s 6 6 (a, 0 p 0 a, n ) = 0 

~ r q 4mf4-2u-r 2u-q 4nrK3-2u-r 2u+l-q 

n=0 'l 'l n 

u=0 



(3.13) 



r ,q=0 

In equation 3.9 when the exponent 2(u+v) = 4m+2 (m=0 , 1 , . . . , (N+T- l)/2 ) , 
the summation constraint becomes eq. 3.12; and when the exponent 2(u+v) = 
4mf4 (m=0 , 1 , . . . , (WT-2)/2) , the summation constraint becomes eq. 3.13 
since the summation must be non-negative when evaluated at s=jo). 

Collecting terms and simplifying, eqs. 3.12 and 3.13 become 



E „ 6 r\ ( “r^e+l ' V 2 V 5 ° 

u,r,q,m,=0 


(3.14) 


S 66 (a 3 -a ,3 , •, ) = 0 

n r q ye 7-I p+1 

u,r,q > m=0 ^ c 


(3.15) 
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\ 




where 



u=0 , 1 , . . . ,N+T/2 
r-0,1,.. . ,T 



u=0 , 1 , . . . ,(NfT-l)/2 




q 0,1,... ,T 




m=0 , 1 , . . . ,NfT-l/2 



> • • • > 



(N*fT-2)/2 



(3.14) 



(3.15) 



and 



7 = 4nrf4-2u-4 0 = 2u-q 



N = Order of original impedance. T = order of augmenting polynomial. 

The final value of the indices is chosen so that if the value is non- 
integer the rule is to round down to the next lower integer. 

In the summations above, given values for a and 3 , 6 f s must be chosen 
so that these inequalities are satisfied. 

Expanding the generalized constraints for a first-order augmenting 
polynomial (i.e., T = 1; r,q = 0,1) 



A6 Q 2 + + C6 1 2 = 0 




(3.16) 



where 



~ Ct '4nrM-2u^2u+l ^4m+2-2u^2u 



~ a 4m-2u^2u+l a 4m+2-2u^2u-l 



~ a 4m-2u^2u a 4nrhl-2u^2u - 1 



(3.17) 



— a 4nrf4-2u^2u a 4m+3-2u^2u+l 



“ a 4irrf4-2u^2u - 1 a 4m+2-2u^2u+l 



— a 4nrf3-2u^2u - 1 a 4m+2-2u^2u 
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a l^l " ^2^0 " a 0^2 






a o 


A = 


a 3 P 3 + a 1^5 + a 3^1 " a 2 P 4 ' a 4 P 2 


' a 0 S 6 ' a 6 S 0 


= 


a l 




• 






• 



(3.18) 



B - 0 (3.19) 





a o p o 




’ c o 


c = 


a 4^0 + a 0^4 + a 2^2 " S 13 ! ’ a l^3 




C 1 



(3.20) 



D = 




(3.20a) 



E = 0 



F = 




(3.20b) 



(3.20c) 



and so for a first-order augmenting polynomial the constraints are 

( 3 . 21 ) 



Vo 2 + c k 6 i 2 a 0 



C k+1 6 0 + a k 6 l 0 



23 



where k = 0, 1 , . . . , (N- l)/2 for the a's and k = 0.1,..., N/2 for 
the c f s; N is the order of the original driving point impedance. 

Performing the expansions above for higher-order augmenting poly- 



nomials provides 


the 


results 


displayed 


in 


Tabli 


The equations of 


Table II can 


be put into i 


3.21 becomes 


















6 0 


•i] 


| 


a k 


0 




5 0 








1 


0 


c k 




6 1 




6 o 


6 i 




C k+1 


0 




6 0 




• 






0 


a k 




6 1 



= 0 



= 0 



(3.22) 



In general it can be shown that 

6 fc 6 > 0 

- 

6 fc 6 > 0 



(3.22a) 



where superscript t indicates a transpose, k = 0,1,..., N/2 and 
j = T+l for the first-order example (T is the order of the augmenting 
polynomial) . 



-c. 



k-1 



0 



k-1 



- a u 0 



'k-i 



'k-1 



0 a 



k-1 



• a k ° 



'k-1 



C k-1 ° 
0 a 



k-2 



a k-2 ° 



'k-1 



-a 



k-2 



'k-2 



k 



(3.23) 
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'k+1 



-a. 



0 -a. 



a. 0 

k 

0 c. 



• c k 0 



0 c, 0 

k 

- c k 0 Vi 

0 -Vi 0 

Vi 0 - C k-1 



k-1 



0 c 



a k-l 0 



k-1 0 



■ C k-1 0 a k-2 



= RV (3.23) 

“J 



k k 

Note that and are dimensioned (T+l) x (T+l) . 

To illustrate the process so far developed consider the dp impedance 
used for the example in chapter V. 



Z (s ) = 



s + l/2s + 1 



s + s + 1 

which requires a third-order augmenting polynomial. 

T = 3 N = 2 



a o = 1 



e 0 = 1 



a L = 1/2 



3l = l 



a 2 = l 

P 2 = l 



j = T + 1 = 4 
k = 0 for a f s 
k = 0, 1 for c 1 s 

a 0 = a l^l " a 2^0 " a 0^2 = _3 ^ 2 



= 0 



C 0 = 



C 1 = 



a o^o 

° 2 e 2 



= 1 
= 1 
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-3/2 


0 


-1 


0 




0 


0 


-1 


0 


0 


1 


0 


3/2 


1 


0 


1 


0 


0 


-1 


0 


0 


0 


^4 = 


-1 


0 


-3/2 


0 


0 


3/2 


0 


0 




0 


0 


0 


1 


1 


0 


3/2 


" 

0 




1 


0 


0 


0 


0 


-3/2 


0 


-1 


1 


0 


0 


0 


-1 


3/2 


0 


1 


0 


^4 = 


0 


0 


1 


0 


0 


-1 


0 


0 




0 


-1 


0 


-3/2 













_ 






— 





It remains to solve for the 6 coefficients according to eq. (3.22) 
where 




B. LINEARIZATION AND SOLUTION OF CONSTRAINTS 3,11,12 

In this section the non-linear inequalities developed in the pre- 
ceding section are linearized using a set of congruent transformations 
and solved. Inspecting matrices (J and R, the pattern of their formation 

becomes apparent. Note that they are symmetric for all orders. Evalu- 
1c 1c 

ating £ and FL for fixed j=T+l over the indices k=0, 1, . . . ,N/2 produces 
a set of constraint inequalities as follow 
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(3.24) 



or 



6 C 1 5 0 

6 l Rj 6 2 ) 



6 C 2° 6 2 0 



(3.25) 



*' Rf 2 6 2 0 



Inequalities 3.25 represent N conditions involving T+l unknowns in 

6. (1) 

Where 



6 = 



Vi 



As a first step toward the simultaneous solution of the set of 



inequalities above, consider the transforms for eq. 3.24. 

6 = x 
6 = N j x 



(3.26a) 



k k k k 

where M. (N.) is the normalized modal matrix of Q. (R.), i.e., the 
J J J J 

k k 

matrix of normalized eigenvectors of £ (R^). These orthogonal 



( 1 ) 



Since one of the 6. can be chosen arbitrarily, there are actually T 
unknowns in 6 . ^ 
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t 1c k 

transforms convert the quadratic forms ^6 &_ (6. 6_) into a linear 

combination of the squares of the coordinates of the vector x with no 



cross products. Alternately the transform can be viewed as reducing 
k k 

the matrices £ (R ^ ) to diagonal form; 



l.e . , 



6 fc £ k 6 



(3.27) 



under the transform 



6 = x 
“J “ 



(3.27a) 



becomes 



t /w k. _.k _,k t . k.-l Jc _,k , 2 2 , 

x (M.) Q. M. x = x (M.) 0. M. x = 1 x_ + l_x_ + 

“ J J “J “ ~ J 3 ~J _ 0 0 11 



+ 1. ,x. 
J-l J-l 



(3.28) 



k t k -1 

since a property of an orthogonal matrix is (M.) = (M.) . The 

k k 

l's in eq. 3.28 are the eigenvalues of Q , and is the orthogonal 
matrix of normalized eigenvectors which are always available. 



Using the transforms above in eq. 3.25 yields 

t /v-k. t _k _,k > n 

x (M.) Q. M. x = 0 
“ J J "J ~ 



t /M k.t „k „k > . 
x (N.) R.N.x=0 

_ — j -j - 



(3.29) 



Substituting 



D k - (mV 2 k M k 
q-j J J -j 

D k - «£) R k N k 
r~j “J "J “J 



(3.29a) 



where the D's are diagonal matrices, 3.29 becomes 



t k > 

x D. x = 0 

~ q-j ~ 

t k > „ 

x D . x = 0 

~ r-j - 



(3.30) 
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Writing the quadratic forms as the product of a row and column vector 



eq. 3.30 becomes 



where 



(row) z = 0 

<rj 

(row) z = 0 

r~J 



z = 



X j 2 -1 



(3.31) 



(3.31a) 



1c 1c 1c 1c 

and D. (row) ( D. (row)) indicate D. ( D.) written as a row. There 

q~j q-j r-j 

are N constraints on T+l variables in £. 

The solution to this linearized set of inequalities can be quite 

2 

elegantly found using Danzig*s simplex method modified to permit 
optimization according to the criteria that as many of the individual 
inequalities as possible set to zero. The simplex method initially 
determines if the set of inequalities is consistent. For the purpose 
of this thesis it is recommended that one start with T = 1. If the 
inequalities are inconsistent simply increase the order of the aug- 
menting polynomial by one and proceed until a feasible set is generated. 

1 

According to Guellemin , if the even part of Z(s) has no jco-axis zeros a 
finite feasible set can always be found. 

A rigorous proof of this assertion is not available in the literature. 
Next the inverse transformation of eq. 3.26a must be accomplished in such 
a way that the resultant 6's are equal (i.e., 6 q=6q, 6^=6^, etc -)« 
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It should be noted that since and are different, different 

-J -J 

6's can be expected for the same x's. In Eq . 3.27 pos tmultiplying the 
k -k k k 

(Hj) b y ^ ( n 2 )? which are arbitrary diagonal matrices, doesn't 

alter the congruency of the transform (i.e., the quadratic form in 6_ is 

still transformed to a linear combination of the squares of the co- 

k k 

ordinates of the vector x) . Since M. (N.) is the normalized modal 

-J -J 

matrix, pos tmultiplying by an arbitrary diagonal matrix is equivalent 
to scaling the modal matrix. 

Therefore, inverse transforming eq. 3.26a in such a way that the 6 f s 

k k 

are equal is equivalent to requiring that diagonal matrices W and W 

m— n— 

be found such that 



-.k TT k __k TT k 
M. W x = N. W x 

- j nr- - -J rr- - 



(3.32) 



where eq. 3.32 is developed by setting the two equations 3.26a, 

postmultiplied by the required diagonal matrices, equal to each other. 

k k Jc k 

For the purposes of this thesis choose = 1 1^ and = p I_, 

k k 

where 1 and p are scalars and 1^ is the identity matrix. With this 
substitution, eq. 3.32 becomes 



,k . k k ..k 

1 M. x = p N. x 
“J ~ -J ~ 



(3.33) 



or 



M* * - q k>k 
“"j * ^ 



X 



(3.33a) 



where 



k,k k k 

q = i p 



,k.t 



Rewriting 3.33a by premultiplying both sides by (N^) yields: 

/vrkvt w k k.k 

(N.) M. x = q ’ x 
“J ~ 



(3.33b) 
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Expanding eq. 3.33b over the index k,(j = T+l) yields 



. O.t 0 0,0 . 0. t 1 0,1 

(Nj) x = q x (KL) x = q x 

,.,Lt 0 1,0 /M l.t 1,1 

(N.) M . x = q x (N.) M . x = q x 

-J ~J ~ - ~J ~J “ 



,„0 N t N/2 

(N.) M. x 
”J “J “ 

, M Ut N/2 
(N.) M. x 
"J "J ~ 



0 ,N/2 

q x 

1 ,N/2 

q x 



(3-34) 



,J/2,tC N/2 ,0 . N/2. t 1 N/2,1 

(NV ) M. x = q ’ x (N. ) M.x = q ’ x 

-J “J “ “ “J “J" 



- q N/2 ' N/2 £ 



Equations 3.34 are the familiar eigenvalue problem in which solution 

k k k t k 

constants q 5 (the eigenvalues) can always be found for (NL) Mj non- 
singular. Since equations 3.34 can always be solved for the eigenvalues, 
eq. 3.26a can always be inverse transformed such that the 6*s are equal; 
and, consequently, a vector 6^ can always be found which satisfies each 
of the constraints in eq. 3.25 simultaneously once a consistent set has 
been generated and solved for the vector x. 

Fortunately it is not necessary to solve for 6^ explicitly. It is 
sufficient to realize that evaluating eq. 3.30 at the determined solution 
vector x yields numerical values: 



t _k . 

x D . x = k. 

- q-j - k 

t _k , 

x j x = k^ 



(3.35) 



t k 

If a particular scalar k^ is zero, then the corresponding (£.6^ 
t k 

(6 R,6) is also zero; and if a particular k is positive, (the only 
J k 

t k t k 

other possible condition) then the corresponding _6 (J 6. (6. R are also 
positive although of different magnitude. This can be demonstrated as 
follows : 



Let 



* T7 k 

6 = M. W x 
- -J nr - 
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which is 



for 



6=1 M. x 
-j - 



= i k i 



(3.35a) 



Substituting eq. 3.35a in 3.24 yields: 

t , k.t ,,k 2 k k 1 . 1 k. 2 2 /ik.2 2 

x (Mj) (1 ) 3 M x = 1 0 (1 ) x Q + lj_ (1 ) x L + . 



,.k.2 ri 2 2 

= (1 ) [l 0 x 0 + l lXl + 



+ ViYi 1 



, ., k % 2 2 
• + Vi <l 5 x j-i 

(3.35b) 



Comparing eq. 3.35b with 3.28, it is evident that eq. 3.35b is eq. 3.28 
multiplied by a positive constant which establishes the assertion above. 

Recalling that the quadratic forms are actually the numerator co- 
efficients of the real part of the augmented dp impedance, by solving 
eq. 3.31 and evaluating 3.30 the corresponding inequalities 3.24 are 
revealed as either zero or positive. It will now be demonstrated that 
this information is sufficient to determine the optimal augmenting and 

augmented polynomial to insure that the real part of Z (s) is positive 

a 

term by term. 



C. OPTIMAL AUGMENTING AND AUGMENTED POLYNOMIALS 

In this section knowledge of the form of the augmented driving-point 
impedance (i.e., which coefficients are positive and which zero) will 
lead to the development of the optimal augmenting polynomial and con- 
sequently the optimal augmented driving-point impedance in a straight- 
forward way. 

Consider an impedance function 



Z(s) = 



m 2 + n 2 



(3.36) 
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and an augmenting polynomial 



P(s) = m 1 ' + n x ' . 



Then 



P( s ) ®i" + “i" 

Z (s) = Z(s) “t = ~ ,, - V - „ 

a v ' v ' P(s) m 2 " + n 2 " 

is the augmented polynomial. 

m "m 1 - n "n " 

ReZ ( joo) = s j s=joo 

(m o"> - (n ") 



where 



•• 

m l 


= 


m^m^ 1 


+ 


n l n l 


m 2* 


= 


m 2 m i 


+ 


Vl 


n i' 


= 


m l n i 


+ 


n l m l 


n 2 " 


= 


m 2 n i 


+ 


n 2 m l 



(3.36a) 



(3.37) 



(3.38) 



(3.39) 



Substituting eq. 

ReZ (jco) = 
a 



3.39 in 3.38 yields 

(n^' + n^) (m 1 1 - n^) (m^ - n^) 
(m^ + n 1 ') (m 1 ' - n^) (m 2 2 - n 2 2 ) 



S=jOO 

(3.40) 



ReZ (jco) = Re (P(ico) ) ReZ (joo) 
P(jco) 



(3.41) 



Num(ReZ a (jco)) ^ R e p(i 0L) ) Num(ReZ (jco) ) . ... 

Den (ReZ (joo)) Den (ReZ (joo)) ^ ) 

a a 

where Num(ReZ (jco)) and Den(ReZ (jco)) represent numerator and 
a a 

denominator terms respectively. Therefore, 

Num(ReZ (joo)) = ReP(joo) Num(ReZ (joo) ) (3.43) 

a 

From eq. 3.24 we know which coefficients of Num(ReZ (jco)) are zero 

a 

and which are positive, and we also know Num(ReZ (jco) ) , the original 
unaugmented dp impedance. A simple-minded but none the less effective 
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way to find ReP(jco) and the unknown coefficients of Num(ReZ (jco)) is 

a 

to proceed by long division. The scheme is best illustrated by a 
numerical example. 

Given Num(ReZ ( jco)) = co** - 2co^ + 3co^ - 4oo^ + 5 , developing the form 

of Num(ReZ (jco)) using the Danzig method outlined above to solve eq. 
a 

3.31 for z and substituting into eq. 3.30, the following information 

is obtained concerning Num(ReZ (jco)) 

a 

Num (ReZ^ ( jco) ) = (Xi^+0a3^^+0a)^+0aj^+0aj^+c^aj^+c^ , 

where c.. > 0 and c > 0 . 

-> o 

Dividing Num(ReZ ( jo)) ) into Num(ReZ (jco)) 

a 

1 +2 +1 

1 -2 +3 -4+5 | 1 +0 +0 +0 +0 +c +c 

1 -2 +3 -4 +5 * 

+2 -3 +4 -5 +C 5 

+2 -4 +6 -8 +10 

+1 -2 +3 +(c,.-10) +c, 

+1 -2 +3 -4 * +5 b 

(c 5 -6) (c 6 -5) 

and setting = 6 and = 5 the division terminates. The numerator 

^ 2 12 2 

of the ReP(jco) = co + 2co +1 and the numerator of ReZ (jco) = co + 6 co +5. 

a 

A feature of this method that can be useful is that the augmenting poly- 
nomial ReP(jco) always has non-negative coefficients. It may even be 
Hurwitz . In the Miyata method this added constraint is not needed; but 
in transfer function synthesis, where augmentation may be called for, 
the augmenting polynomial is additionally required to be Hurwitz. 

It turns out that if Z(s) is of second order or less, it is unnecessary 

4 

to have exact knowledge of the form Num(ReZ (jco)) since Poincare and Lewis 

a 

have demonstrated that for Z(s) quadratic an augmented polynomial having 
but three non-zero positive terms can invariably be generated. For Z(s) 
of greater than second order, in general, exact knowledge of the form of 
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Num(ReZ (jco)) is required and the preliminary steps which yield this 
a 

form cannot be bypassed. 

The problem is now essentially solved, and a brief outline of the 
computational steps is in order. 

First generate a consistent set of non-linear inequalities (3.24) 

and then employ the transform 3.26a to set up and solve the resultant 

linear set of inequalities 3.31 by the Danzig simplex method which 

insures that as many as possible of the inequalities are zero, thereby 

minimizing the number of elements in the subsequent network synthesis. 

Armed with the solution vector x> determine the form of the solved 

inequalities 3.35 that yields which coefficients of the numerator of 

the real part of the augmented dp impedance are zero and which are 

positive. Divide the given Num(ReZ (jo)) ) into Num(ReZ ( jcu) ) and solve 

a 

for the unknown positive coefficients. The process is now complete and 
the numerator of the real part is known (evaluating the denominator is 
straightforward), and one can proceed to the synthesis part of the 
problem with a polynomial that is positive term by term. 

At this point it seems wise to demonstrate the solution algorithm 
and digital computer program developed to augment a moderately involved 
driving-point impedance . 
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IV. DIGITAL COMPUTER PROGRAM 



From consideration of the computational scheme outlined, it is 
obvious that the technique is admirably suited to implementation on a 
digital computer. This is espacially true since solving a set of linear 
inequalities under the constraint that the solution variables be greater 
than zero is very similar to the basic linear programming problem which 
is the subject of the Mathematical Programming System"* (MPS/360). 

Therefore, a highly developed and efficient system is readily 
available to solve the linearized inequalities, and the programming 
problem becomes one of simply generating the inequalities from the 
original driving-point impedance and inputing the set into the MPS 
system. 

The digital program required for systematic solution is straight- 
forward and, aside from typical programming difficulties, simple. The 
entire solution process has not been consecutively programmed due to 
interface restrictions between OS/360 and MPS/360. The program follows 
in block diagram and Fortran format. 
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SUBROUTINE GROAN 



Calculate the element values for the 



V 

SUBROUTINE 



V 

SUBROUTINE 



V 

SUBROUTINE 



V 

MPS/360 



<- 



matrix of nonlinear inequalities. 



GRUNT 



Form the matrix of inequalities. 



GROWL a Convert the matrix of nonlinear 

inequalities to a matrix of linear 
inequalities . 



PUNCH 



Prepare an input data deck for the MPS/360 
program which solves the linear inequalities 
and consequently the form of the original 
nonlinear inequalities . 



n=n+l 

-> 

Output the required order of augmentation 
and form of the numerator of the real part 
of the driving-point impedance. 
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// EXEC FORTCLGPfPEGIONo GO=l OOK, T I MF» GQ= 5 
//FORToSYSIN DD * 

DIMENSION ALFA( 10) ,BETA( 10 ) , A ( 10 ) , C ( 10 ) , 0 ( 10 , 1 0 , 1 0 ) 
DIMENSION R(lOtlOtlO) 

DIMENSION QR(10,10) , XX (10,36 ) ,YY( 10,36) 

1 FORMAT (1 10) 

2 FORMAT ( 8F10o2 ) 

READ (5,1) N 
READ (5,1) N 
M3= (M+2)/2 
M4=2^M3 
M5=4*M3 
MM=M+1 
N1=N+1 
N2=N1+M 

READ ( 5,2) ( ALFA ( I ) ,1=1, MM) 

READ (5,2) (BETA(I ),I=1,MM) 

CALL GROAN ( M , N2 , M3 , M5 , ALF A , BETA , A , C ) 

D03 N9=l ,N1 

CALL GRUNT ( M3 , N9 , A , C , Q , R ) 

CALL GROWL ( Q, R , M3 , M4 , N9 , QR , XX , YY ) 

CALL PUNCHY (N2,N9,M4»QR) 

3 CONTINUE 
STOP 
END 



SUBROUTINE GROAN ( M,N2 , M3, M5 , ALFA, BETA, A, C ) 
DIMENSION ALFA(IO) , BETA ( 10 ) , A( 10 ) , C ( 10 ) 

MM=M+2 

DO 1 I =MM , M5 
ALFA(I )=0 o 0 

1 BETA ( I ) =OoO 
D02 Ml = 1 ,M3 

A ( Ml )=OoO 

C( Ml )=0o0 
M2=2*M1 
DO 3 N1=1,M2 
I=2*N1 
J = I - 1 
K=4* Ml- 1 
L = K+ 1 
11=1-1 
JJ=J-1 
KK=K-1 
LL=L-1 
I F ( K ) 4,4,5 

4 K = MM 

5 IF(KK) 6,6,7 

6 KK=MM 

7 IF (LL) 8,8,9 

8 LL=MM 

9 IF (JJ) 10,10,11 

10 JJ=MM 

11 A( Ml ) = ALFA ( I)X'BETA(K)-ALFA(J) 5 ) t BETA(L)+A(Ml) 

C ( Ml )= ALFA(II )*BETA(KK)-ALFA( J J )*BETA( LL ) +C ( Ml ) 
3 CONTINUE 

2 CONTINUE 
M6=M3+1 
D012 I = M 6 , N 2 
A ( I ) =0 o 0 

12 C ( I ) =Co 0 
RETURN 
END 



SUBROUTINE GRUNT ( M3 , N9 , A , C , Q, R ) 

DIMENSION A( 10 ) ,C( 10) ,Q( 10,10, 10) ,R( 10,10,10) 
DO 1 1=1, M3 
DO 2 J = 1 ,N9 
DO 3 K=1,N9 
Q ( I , J, K ) =0 o C 
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3 R ( I » J, K ) =0 o C 
2 CONTINUE 

1 CONTINUE 
DO 4 1=1, M3 
L = I 
N = 2* I 

DO 5 J=1 ,N , 2 
DO 6 K= J , N9 , 4 
Q( I, J,K)=A( L) 

Q ( I , J+l , K+ 1 ) =C ( L ) 

Q ( I , Jt K+2 ) =-C ( L ) 

IF (J-l) 7,7,11 

7 Q ( I , J + l , K+3 ) =-A ( L ) 

GO TO 12 

11 Q( I , J+l , K+ 3 ) =-A ( l-l ) 

12 IF(IoEQol) GO TO 4 
IF(L-l) 4,5,6 

6 L=L- 1 
5 CONTINUE 

4 CONTINUE 

DO 15 1=1, M3 

L= I 

N=2*I 

DO 16 J= 1 , N , 2 
DO 17 K= J , N9 , 4 
R ( I , J,K)=C(L+1) 
R(I,J+1»K+1)=A(L) 
R(I,J+2,K+2)=C(L) 

R ( I , J, K+2) =-A ( L) 

R ( I , J+l , K+3 ) =-C ( L ) 

I F ( I o EQo 1 ) GO TO 15 
IF (L-l) 15,16,17 
17 L=L-1 
16 CONTINUE 
15 CONTINUE 
DO 8 1=1, M3 
DO 9 J=1,N9 
M = J+ 1 

DO 10 K=M,N9 
Q ( I , K, J ) =Q ( I , J , K ) 

10 R ( I » K, J ) =R ( I , J , K ) 

9 CONTINUE 

8 CONTINUE 

DO 20 1=1, M3 
DO 30 J=1,N9 

30 WRITE<6,25) (0(1, J,K) ,K=1,N9) 
DO 40 J=1,N9 

40 WRITE(6,25 ) <R( I, J,K) ,K=1,N9) 
20 CONTINUE 
25 FORMAT (//, 6F15o5) 

RETURN 

END 



SUBROUTINE GROWL ( Q , R , M3 , M4 , N 9, QR , XX , YY ) 
DIMENSION Q (10, 10, 10 > ,R( 10,10,10) 
DIMENSION QR(10,10) , XX (10,40 ),YY( 10,M)) 
DIMENSION Q C (10) , RR ( 10 ) , X ( 40 ) , Y ( 40 ) 

DO 12 1=1, M3 
11 = 1 

DO 9 K=1,N9 
DO 10 J=1 , K 
QQ ( I I ) =0 ( I , J,K) 

RR ( I I ) =R ( I , J,K) 

10 11=11+1 
9 CONTINUE 

CALL EIGEX( CO, X,N9,0) 

CALL EIGEX(RR,Y,N9,0) 

L = 1 

I I 1 = 2*1 -1 
I 1 1 1 = 1 1 1+1 
DO 15 J = 1 1 NS 
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QR ( I II , J)=QQ(L) 

QRUIII, J)=RR(L) 

15 L= J+L+l 
N10=N9**2 

DO 16 J=l, N10 
XX( I » J)=X( J ) 

16 YY ( I , J ) = Y ( J) 

12 CONTINUE 

DO 40 1=1, M4 

40 WRITE(6,30) (QR(I,J) ,J=1,N9) 
DO 45 1=1, M3 

WRITE (6, 35 ) ( XX ( I , J ) , J = 1,N10 ) 
45 WRITE<6,35) (YY( I,J) ,J=1,N10) 
30 FORMAT (/ /»6F15o5) 

35 FORMAT (//,6F15o5) 

RETURN 

END 



SUBROUTINE PUNCHY ( N2 , N9 , M4 , QR ) 

DIMENSION QR(10, 10) ,QR1( 10,10) 

40 FORMAT (• NAME’ ,10X, • TESTPROB’ ) 

41 FORMAT ( » ROWS' ) 

42 FORMAT ( 2X , ' N COST • ) 

43 FORMAT ( 2 X, 'G ROW', ID 

44 FORMAT ( 'COLUMNS' ) 

45 FORMAT (4X, 'C0L1' ,6X, 'COST' ,6X, 'OoO' ) 

46 FORMAT (4X, 'C0L1' ,6X, 'ROW* , II ,6X , F12o 5 ) 

47 FORMAT ( ' RHS ' ) 

48 FORMAT ( 4X ,' RHS ', 7X ,' COST ', 6X ,' 0 o 0 • ) 

49 FORM AT (4X, 'RHS' ,7X, 'ROW' ,Il,6X,F12o5) 

50 FORMAT ( • ENDATA* ) 

60 FORMAT (4X, • COL' , II ,6X, 'COST* ,6X, 'OoO' ) 

61 FORMAT (4X, 'COL', I1,6X, 'ROW' , I 1 , 6X, F12* 5) 

62 FORMAT (4X, 'RHS',7X, • ROW • , 1 1 , 6X , F12o 5 ) 

I F ( N9o EQo 1 ) GO TO 70 

IF ( N9o GT© 2 ) GO TO 30 
WRITE! 7,40 ) 

WRITE ( 7 , 41 ) 

WRITE(7,42 ) 

WRITE! 7,43 ) (I ,I = 1,M4) 

WRITE! 7,44) 

WRITE(7,45 ) 

WRITE(7,46) (I,QR(I,1),I=1,M4) 

WRITE! 7, 47 ) 

WRITE(7,48 ) 

DO 1 1=1, M4 

1 QR1 ( I , 2 ) =-QR ( 1,2) 

WRITE! 7,49) ( I,QR1 ( I ,2) , 1=1 ,M4) 

WRITE! 7, 50 ) 

GO TC 70 
30 WRITE ( 7 , 40 ) 

WRIT E( 7, 41 ) 

WRITE(7,^2 ) 

WRITE! 7,43) ( I ,I=1,M4) 

WR IT E ( 7 , 44 ) 

N8=N9-1 

WR IT E( 7 , 60 ) ( 1=1, N8 ) 

D02 J=1,N8 

2 WRITE! 7,61) (J,I,QR(I,J) ,I=1,M4) 

WR IT E ( 7 , 47 ) 

WR I T E ( 7 » 48 ) 

D03 1=1, M4 

3 QR1 ( I , N9 ) = -QR ( I , N9 ) 

WRITE(7 ,62 ) (I,QR1( I,N9),I=1,M4) 

WRITE! 7, 50) 

70 RETURN 
END 
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SUBROUTINE E I GEX ( A , R, N , MV ) 

DIMENSION A <1 > ,R <1 ) 

5 RANGE=lo 0E-6 

IF(MV-l) 10,25,10 
10 IQ=-N 

DO 20 J=1,N 
I G= I Q+N 
DO 20 1=1, N 
I J= I Q+ I 
R ( I J ) = 0o C 
I F ( I - J * 20,15,20 
15 R< IJ ) = lo0 
20 CONTINUE 
25 AN0RM=0o 0 
DO 35 1=1, N 
DO 35 J = I , N 
IF(I-J) 30,35,30 
30 I A= I +( J* J- J )/2 

ANORM=ANORM+A( IA)*A( I A) 

35 CONT INUE 

IF(ANORM) 165,165,40 
40 A NORM= lo 41 ASSORT ( ANORM ) 

ANRMX= AN OR MAR ANGE/ FLOAT ( N ) 

I ND=0 
THR= ANORM 

45 THR=THR/FLOAT(N) 

50 L=1 
55 M = L + 1 
60 MQ= ( M*M-M ) /2 
LQ= ( L*L-L ) /2 
LM=L +MQ 

62 I F ( ABS( A( LMM-THR ) 130, 65,65 
65 I ND=1 
LL=L +LQ 
M M=M+MQ 

X=0o5*(A(LL)-A(MM) ) 

68 Y=-A(LM)/ SQRT(A(LM)AA(LM)+XAX) 

I F ( X ) 70,75,75 
70 Y=- Y 

75 S INX=Y / SQRT(2oO*( loO+< SQRT ( lo 0-Y*Y ) ) )) 
SINX2=SINX*SINX 
78 C0SX= SQRT ( lo 0-S I NX2 ) 

COSX2=COSX*COSX 
SINCS =S I NX+COSX 
ILQ=N*( L-l ) 

I MQ=N# ( M-l ) 

DO 125 1=1, N 

i q= ( i*i-n / 2 

IF(I-L) 80,115,80 
80 IF(I-M) 85,115,90 
85 I M= I +MQ 
GO TO 95 
90 I M = M + I Q 

95 IF(I-L) 100,105,105 
100 I L= I +L Q 
GO TO 110 
105 IL=t+IQ 

110 X = A< IL)*COSX-A(IM)*SINX 

A( I M ) =A ( IL ) * S I NX + A ( I M ) ACOSX 
A ( IL ) = X 

115 IF(MV-l) 120,125,120 
120 I LR= ILQ+ 1 
I MR= IMQ+I 

X = R< ILR)*CCSX-R( IMR)*SINX 
R( I MR) =R( I LR)*SINX + R< IMR)*COSX 
R ( I LR I =X 
125 CONTINUE 

X=2oO*A(LM)*SINCS 
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Y=A( LL)*C0SX2+A(MM )*SINX2-X 
X=A ( LL)*SINX2+A(MM)*COSX2+X 

A( LM ) = ( A (LL )-A(MM) )*SINCS+A( LM ) *( COS X2- S I NX2 ) 
A ( LL ) =Y 
A ( MM ) = X 

130 IF(M-N) 135,140,135 
135 M=M+ 1 

GO TO 60 

140 I F ( L-( N-l ) ) 145,150,145 
145 L=L+1 

GO TO 55 

150 IF(IND-l) 160,155,160 
155 I ND=0 

GO TO 50 

160 I F ( T HR-ANRMX ) 165, 165,45 
165 RETURN 
END 

//GCoSYSIN DD * 
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//JOBLIB DO DSNAME=SYSl o MPS360LP,DISP=(SHR, PASS) 
// EXEC L INPRGG 

//MPSloSYSIN DD * 



INFEAS 


PROGRAM 
I NIT I AL Z 

MOVE ( XPBNAME, • PBFILE' ) 
MOVE ( XDA7A , 'TESTPROB' ) 
MOVE ( XO B J , 'COST' ) 

MOV E ( XRHS » ' RHS ' ) 
MVADR(XOONFS, INFEASI 
MVAOR(XDOFEAS* FEAS ) 
GCTO( AARRRG) 

MOVE ( XPBNAME, 'PBFILE' ) 
MOVE ( XDATA, 'TESTPROB' ) 
CONVERT ( • SUMMARY* ) 
GOTO( AAARRG) 


AARRRG 

AAARRG 


CONVERT ( 'SUMMARY* ) 

BCOOUT 

SETUP! ' MAX' ) 

PRI MAL 


FEAS 


SOLUTION 

EXIT 

PEND 



//MPS2o SYS IN OD * 
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V. SAMPLE PROBLEM AND SOLUTION 



zoo - 4 + 1/2s + 1 

s + s + 1 

From computer program results 

NRez^Cjco) = cu 10 +0oo^ +0cjd^ + Oco^ +c^co 2 +c^ 

NRez(jco) = co^ - 3/2co 2 + 1 

DReZ(jco) = co^ -co 2 + 1 
Dividing NReZ (jco) by NReZ(jco) gives: 

3 



1 +3/2 +5/4 +3/8 

1 -3/2 + 1 1 0 0 0 +c, +c 

1 -3/2 1 4 * 

3/2 -1 

3/2 -9/4 +3/2 

+5/4 -3/2 

+5/4 -15/8 +5/4 

+3/8 -5/4 

+3/8 -9/16 +3/8 

-11/16 -3/8 

Choosing and c,_ equal to 11/16 and 3/8 respectively, the division 
terminates and the real-part numerator of the augmenting polynomial is 
seen to be 



co ^+3 / 2o/*+5 / 4 co 2 +3 / 8 . 



The real-part numerator of the augmented polynomial has three non- 
zero terms and reads 



10 2 
£0 +ll/16a) +3/8. 



The denominator is 



DReZ(jco) x ReP(jco) 



which equals 

to^+1/ 2co^+3 / 40)6+5/ 8co^+7 / 8co 2 +3 / 8 . 
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Therefore 



ReZ a (jcu) 



1 0 ? 

to +ll/16cu +3/8 

0) 1 °+ 1 / 2cd 8 +3 /4u 6 + 5 / 8co 4 +7 / 8 cd 2 +3 / 8 



which is positive for all co, as required. 
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TABLE I 6 

COMPARISON OF THE NUMBER OF CIRCUIT ELEMENTS REQUIRED FOR SINGLE-N 





AND 1/2 -N 


SPLIT 


MIYATA AND MODIFIED 


BOTT-DUFFIN METHODS 




Single-n 

Split 






Half-n 

Split 


Modified 
Bott-Duf f in 
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2 
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2 
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2 


2 


4 


3 


4 


3 


5 


3 


3 


9 


4 


6 


4 


7 


5 


4 


16 


5 


10 


5 


15 


7 


5 


25 


6 


14 


6 


19 


11 


6 


36 


7 


17 


7 


35 


15 


8 


64 


9 


25 


9 


75 


31 


10 


100 


11 


35 


11 


155 


63 


12 


144 


13 


44 


13 


315 


127 


14 


196 


15 


52 


15 


635 


255 


16 


256 


17 


62 


17 


1275 


511 


18 


324 


19 


74 


19 


2555 


1023 


20 


400 


21 


86 


21 


5115 


2047 


30 


900 


31 


139 


31 


163,835 


65,535 



= number of reactive components 
n^ = number of resistive components 
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GENERAL POLYNOMIAL EXPANSIONS OF FIRST THROUGH FIFTH ORDER AUGMENTATIONS 
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